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Abstract. We report the results of a recent (July 2004) XMM-Newton Target of Opportunity observation of the Anomalous 
X-ray pulsar IE 1048.1-5937 together with a detailed re-analysis of previous observations carried out in 2000 and 2003. In 
July 2004 the source had a 2-10 keV flux of 6.2x 10“^^ erg cm“^ s“^ and a pulsed fraction Pf= 0.68. The comparison of the 
three data sets shows the presence of an anti-correlation between flux and pulsed fraction, implying that previous estimates of 
the source energetics based on the assumption of a large and constant pulsed fraction might be significantly underestimated. 

The source spectrum is well described by a power law plus blackbody model (kT~0.63 keV, photon index r~2.7-3.5) or, 
alternatively, by the sum of two blackbodies of which the hotter is Comptonized by relativistic electrons. In this case the 
temperatures are kTi ~0.2-0.3 keV and kT 2 ~0.4-0.5 keV and the emitting area of the cooler component is consistent with 
the whole neutron star surface. The long term luminosity variation of a factor ^2 is accompanied by relatively small variations 
in the spectral shape. Phase resolved spectroscopy indicates a harder spectrum in correspondence of the pulse maximum. No 
spectral features have been detected with Ao limits on the equivalent width in the range ~ 10-220 e’V, depending on line energy 
and width. 
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1. Introduction 

The Anomalous X-ray Pulsars (AXPs, Mereghetti & Stella 
1995; van Paradijs, Taam & van den Heuvel 1995) were origi¬ 
nally considered a subclass of the accreting X-ray pulsars char¬ 
acterized by spin period of a few seconds, stable spin-down, 
soft X-ray spectrum and luminosity larger than the rotational 
energy loss. The absence of bright optical counterparts and of 
Doppler modulations of their pulse period has then led to the 
conclusion that they are very likely isolated neutron stars (see 
Mereghetti et al. 2002 for a review). Although the possibility 
that they are powered by accretion from a residual disk cannot 
be completely excluded, the model currently considered most 
successful involves highly magnetized neutron stars (“magne- 
tars”), in which the decay of an intense magnetic field (10^"*- 
10^® G) powers the X-ray emission. This model was first pro- 
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posed (Duncan & Thompson 1992) to explain a different class 
of enigmatic objects, the Soft Gamma-ray Repeaters (SGRs) 
and then extended to the AXPs due to their many common 
properties (Thompson & Duncan 1996), including the recent 
observations of bursting activity in a couple of AXPs (Gavriil, 
Kaspi & Woods 2002; Kaspi et al. 2003). 

IE 1048.1-5937 is a key object to understand the con¬ 
nection between AXPs and SGRs and to study the physical 
processes involved in the X-ray emission from these sources. 
Besides being the first AXP from which X-ray bursts were dis¬ 
covered, other properties make IE 1048.1-5937apossible tran¬ 
sition object between AXPs and SGRs: it is one of the AXPs 
with the hardest X-ray spectrum and with the most variable pe¬ 
riod evolution, which are typical characteristics of the two best 
studied SGRs (SGR 1900 h- 14 and SGR 1806-20). 

Elux variations were reported for IE 1048.1-5937 in the 
past, comparing measurements from different satellites, which 
could be affected by large systematic uncertainties (Oostebroek 
et al. 1998). Only recently, the fact that IE 1048.1-5937 is 
a variable X-ray source has been firmly established, when 
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very different source intensities were registered in two XMM- 
Newton and two Chandra observations (Mereghetti et al. 
2004). This variability has been confirmed by the monitoring 
program carried out with Rossi-XTE which showed two ex¬ 
tended flares in its pulsed flux* (Gavriil & Kaspi 2004). Since 
the distance of IE 1048.1-5937 is not known (see e.g. Gaensler 
et al. 2005), we assume a distance of 5 kpc. Following this 
assumption, the 2-10 keV luminosity of IE 1048.1-5937 has 
been observed to vary between ^lO*^^ and ^10^^ erg s“*. 

The XMM-Newton satellite has observed IE 1048.1- 
5937 in three occasions: a short snapshot observation was done 
in December 2000 (Tiengo et al. 2002); a 10 times longer ob¬ 
servation was performed in June 2003, when the source was 
more than a factor 2 brighter (Mereghetti et al. 2004); finally, 
a Target of Opportunity observation was requested in response 
to the detection of a burst from the direction of IE 1048.1- 
5937 on 2004 June 29 (Kaspi et al. 2004). These observations, 
spanning an interval of more than three years and catching the 
source at different intensity states, offer the possibility of a sys¬ 
tematic study of the long term changes in the source proper¬ 
ties based on a homogeneous set of data. Furthermore the high 
statistics obtained in the June 2003 observation allows the most 
detailed phase resolved spectral analysis and search for spectral 
lines ever carried out for this source. 

2. Observations and data analysis 

Table Ogives a log of the three XMM-Newton observations of 
IE 1048.1-5937 (hereafter observations A, B and C). The re¬ 
sults of obs. C are reported here for the first time, together with 
a more detailed analysis of the previous two observations. We 
concentrate mostly on data from the EPIC PN (Striider et al. 
2001) and MOS (Turner et al. 2001) cameras since only in 
obs.B the Reflection Gratings Spectrometer (RGS, den Herder 
et al. 2001) provided a sufficient number of source photons. 
In all the observations the medium thickness optical blocking 
filter was used for the three EPIC cameras, while different ob¬ 
serving modes were used. In particular, during obs. C the PN 
and MOS2 cameras were operated in Timing mode, which has 
a better time resolution (0.03 ms for PN and 1.5 ms for MOS), 
but only mono-dimensional imaging capability, which results 
in a higher background in the source extraction region. 

All the data have been processed with the same version of 
the analysis software (SAS 6.0.0) and using the most recent 
calibration files. Time intervals with high particle background 
were filtered out. This reduces the deadtime-corrected exposure 
time to 16 ks for PN and 15 ks for MOS in obs. C and 40 ks for 
PN and 46 ks for MOS in obs. B. Due to the moderate ampli¬ 
tude of particle background fluctuations and the short exposure 
time, no time filter was applied to the data of obs. A. 

For all the data taken in imaging modes (i.e., PN and MOS 
data of obs. A and B and MOSl data of obs.C), a circle of 30" 
radius centered at the position of IE 1048.1-5937 was used for 

' due to the difficulty of background subtraction and the presence 
of the bright and variable X-ray source Eta Carinae in the field of 
view of this non-imaging satellite, the RXTE data cannot give accurate 
information on the total (pulsed plus unpulsed) flux level. 


the source extraction region for spectral and timing analysis. 
For the Timing mode data the source events were selected from 
a 40" wide strip around the source position. The regions for the 
background subtraction were chosen from the same chip where 
the source was detected. 

No bursts from the source were found by a visual inspec¬ 
tion of the source and background light curves of the three 
observations. This is not surprising because, although EPIC 
has the sensitivity to detect bursts like the three observed by 
RossiXTE up to now, their recurrence time is longer than the 
total exposure time of our observations. 

2.1. Spectral analysis 

For each observation, the 0.6-10 keV spectra from the three 
EPIC cameras were rebinned in order to have at least 30 counts 
per channel and fitted simultaneously using the XSPEC pack¬ 
age vl 1.3.0. To account for the cross-calibration uncertain¬ 
ties between different detectors and operating modes, the rel¬ 
ative normalization was let free to vary. Fixing the PN nor¬ 
malization factor to 1, the values for the MOSl and MOS2 
were 1.02±0.03 and 0.97±0.03 for obs. A, 1.008±0.008 and 
1.043±0.009forobs. B and 1.13±0.02and0.89±0.02forobs. 
C. The larger difference of the last observation values reflects 
the greater uncertainties in the flux calibration of the timing 
mode used in the MOS2 and PN detectors. 

The statistical quality of the obs. B spectrum is so high 
(>130,000 counts in the PN only) that the deviations from 
the best-fit models are dominated by the systematic uncertain¬ 
ties in the instrumental calibration. Therefore we added in all 
the model fits to this data set a 2% systematic error, perfectly 
consistent with the current uncertainty in the EPIC calibration 
(Kirsch et al. 2004). In this way formally acceptable values 
are obtained for the good fits. 

Consistently with previous results on this AXP we found 
that simple models based on a single spectral component are 
not adequate. Therefore we tried different more complex mod¬ 
els, obtaining the best fit spectral parameters summarized in 
Table 2. The “canonical” AXP model consisting of an absorbed 
power law plus blackbody provides a good fit to all the spec¬ 
tra, with absorption N//~ 10^^ cm“^, blackbody temperature 
kTBB~0.6 keV and photon index E in the range ~2.7-3.5. 
The source luminosity was different in the three observations, 
with 2-10 keV observed fluxes of 4.7x10“*^ erg cm“^ s“*, 
10.0x10“*^ erg cm“^ s“*, and 6.2x10“*^ erg cm“^ s“*, re¬ 
spectively in obs. A, B, and C. 

The sum of two blackbody components, successfully used 
in the past for other AXP spectra (e.g., Israel et al. 2001; 
Halpern & Gotthelf 2004), yields acceptable results for obs. 
A, but is clearly rejected by the higher quality spectra of the 
two longer observations. 

In an attempt to physically link the power law and black¬ 
body like components, we also tried a simple Comptonization 
model^. The basic idea is that soft thermal photons, possibly 
produced at the star surface, may be upscattered by relativistic 

^ the applicability of this model to the magnetar case is discussed in 
Sect. 3 
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Table 1. XMM-Newton observations of IE 1048.1-5937 


Observation 

Date 

Duration 

PN Mode'®-’ 
Net exposure 

MOSl Mode^“J 
Net exposure 

MOS2 Mode^“J 
Net exposure 

Period 

Pulsed fraction^*’^ 

A 

2000 Dec 28 

8 ks 

EE 

SW 

SW 

6.45253±0.00009 s 




4.5 ks 

7ks 

7ks 

(91±3)% 

B 

2003 June 16 

50 ks 

EE 

SW 

SW 

6.454894±0.000002 s 




40 ks 

46 ks 

46 ks 

(54.6±0.6)% 

C 

2004 July 08 

30 ks 

Ti 

SW 

Ti 

6.456109±0.000005 s 




16 ks 

15 ks 

15 ks 

(68±1)% 


(a) pp _ p^p pj-ajjjg (time resolution 73 ms); SW = Small Window (time res. 0.3 s); Ti = Timing (time res. 0.03 ms (PN), 1.5 ms (MOS)) 
in the 0.6-10 keV energy range. 


Table 2. Results of phase averaged spectroscopy^®^ 


Model 

Parameter 

A 


C 

PL-tBB 

N/r (lO^^^ cm^'') 

0.95±0.09 

1.08±0.02 



kTss (keV) 

0.63±0.04 

0.627±0.007 

0 f^oTi+o.oos 

U.OZ:>_0.006 


Rbb (km)(") 

0.8±0.1 

1.29±0.03 

1 04+®-^^ 

05 


r 

2.9±0.2 

3.27±0.05 

a 44+0.09 


PL norm^'*) 

3.8±0.3 


9.4±0.3 


Xred (d.O.f.) 

0.963 (255) 

1.046 (519) 

1.041 (432) 

BBl -I-BB2 

Nh (lO^"" cm^'') 

u.jj_0 05 

0.62±0.01 

0.67±0.02 


kTBBi (keV) 

0 47+0-05 

U.H-/_0.06 

0.44±0.01 

0.37±0.02 


Rbbi (km)(") 
kTBB 2 (keV) 

i-Oig;? 

2.8±0.1 

0.86±0.02 

3 

•^•^-0.3 

0 76+0 03 

/0_Q Q2 


Rbb 2 (km)(") 

0.3±0.1 

0.63±0.05 

0.7±0.1 


Xred (d.O.f.) 

1.004 (255) 

1.365 (519) 

1.118 (432) 

CBB 

Nb (lO""^ cm-"") 

0.53±0.04 

0.588±0.008 

0.57±0.02 


kTss (keV) 

0.40±0.02 

0.412±0.004 

0.40T0.01 


Rss (km)("> 

1.7±0.1 

2.75±0.04 

2.3±0.1 



3.8±0.2 

4.40±0.06 

4.4±0.2 


xled (d-o.f.) 

0.994 (256) 

1.401 (520) 

1.273 (433) 

BB-hCBB 

Ns (10'''' cm^'') 

0.8±0.2 

0.82±0.02 

Q 7Q+U.iU 

^•/^_0.05 


kTss (keV) 
Rss (km)("> 
kTcss (keV) 

0 44+0-09 

U.44_o.o5 

0.23±0.02 

8±2 

'J.ZO_o 04 

5-9^14 

0.48t°+ 


Rcss (km)^“^ 

1.5±0.4 

2.3±0.2 

1-612:? 



4 1 +0-7 

+•^-0.4 

4.9±0.1 

5.4T0.5 


Xled (d.O.f.) 

0.964 (254) 

1.026 (518) 

1.042 (431) 


Errors are at the 90% c.l. for a single interesting parameter 
A 2% systematic error has been applied to the model 
Radius at infinity assuming a distance of 5 kpc 

Normalization of the power law component in units of 10“® photons cm“^ s“^ keV~^ at 1 keV 
Comptonization parameter as defined in the text 


e“ (or e^) of small optical depth and mean Lorentz factor ( 7 ). 
An analytical, approximated expression which relates the inci¬ 
dent to the emergent intensity is given in Rybicki & Lightman 
(1979; see also references therein) for monochromatic inci¬ 
dent radiation of energy E', Ie(E) ~ Ii(E')(E/E')^~“. In 
the previous expression a = 1 — In / In A, where 
is the scattering depth in a magnetized medium and A ~ 
4 ( 7^)/3 is the mean energy amplification factor per scat¬ 
tering. The case of a blackbody photon input is easily ob¬ 
tained by convolving the monochromatic expression with the 


initial photon distribution. The photon spectrum is given by 
CE““ dE' E'^+“/[exp (E'/kTBB) - !]■ The model param¬ 
eters are then the blackbody temperature, normalization C and 
the photon index a. We find that the fit with this model alone 
(CBB in Table 2) is not satisfactory, but the addition of a 
second, softer blackbody component, with kT~0.2-0.3 keV, 
results in acceptable fits for all the observations. The hotter 
Comptonized blackbody has a smaller temperature and larger 
emitting area than that found in the blackbody plus power law 
model. 
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Fig. 1. Ratios between the spectra of the three observations (A: top, B: middle, C: bottom) and the (renormalized) best fit model 
of observation B. The colors indicate the three EPIC cameras: PN in black, MOSl in red, MOS2 in green. The data have been 
graphically rebinned to show more clearly the spectral differences. 


Independent of the adopted spectral model, it is clear from 
the best fit parameters reported in Table 2 that the spectrum 
changed between the different observations. This can also be 
seen from FiglU which shows the ratios between the spectra of 
the three observations and the best fit model of obs. B scaled 
only in overall normalization (the normalization factors are re¬ 
ported in the first line of Table 3). The spectrum of obs. C is 
significantly softer than the template spectrum, while that of 
obs. A is slightly harder, although, due to its poorer statistiscs, 
the fit is formally acceptable, as reported in Mereghetti et al. 
(2004). 

For the sake of simplicity we tried to interpret the long term 
spectral variations using the model with less parameters, i.e the 
blackbody plus power law. In order to check whether the long 
term spectral variation could be ascribed to a change in only 
one or two of the four^ relevant parameters we fitted the spectra 
of obs. A and C keeping some of them fixed to the best fit val¬ 
ues of obs. B. The results are summarized in Table 3, where one 
can see that the variations are described in an almost equivalent 
way by a change either in the blackbody temperature (which is 
slightly preferred for obs. A) or in the photon index (which 
gives a better fit to obs. C). If the photon index and tempera¬ 
ture are kept fixed and only their relative normalization or the 
Nh free to vary, a significant improvement in the value is 
obtained for obs. C but not for obs. A. 

The possibility that one of the two spectral components re¬ 
mained constant (both in shape and normalization) throughout 


^ the overall spectral shape, beside depending on Nh, T and kTss, 
is a function of the relative normalization of the two components 


the three observations was also explored, with negative results. 
In fact, due to the flux variation between the two observations, 
fitting the spectrum of obs. A with the parameters of either the 
power law or the blackbody component fixed to those of obs.B 
gives Xred >2 (even if the Nh is left free to vary). 


Table 3. Fits to obs. A and C with a power-law plus blackbody 
model and different sets of fixed parameters. 


Parameter^“^ 

A 

C 

Normalization factor 

0.432±0.006 

0.658±0.004 

xLd (d.o.f.) 

1.081 (261) 

1.341 (438) 

Nif (10^"" cm^^) 

L07±0.02 

L02±0.01 

Normalization factor 

^•^^^-0.009 

0.6321°;°°® 

xLd (d.o.f.) 

1.084 (260) 

1.112 (437) 

Rss (km) 

0.84±0.03 

0.96±0.02 

PL norm 

5.6±0.2 

9.3±0.2 

Xred (d.o.f.) 

1.084 (260) 

1.126 (437) 

kTijij (keV) 

Rss (km) 

0.69±0.02 

0.66l°:“ 

0.602±0.008 

1 07+°°'* 

PL norm 

6.0±0.2 

9.0±0.2 

\red (d.o.f.) 

0.978 (259) 

1.073 (436) 

r 

Rss (km) 

PL norm 

J.OZ_g Qg 

0.72±0.05 

3-0_o.3 

2 /|/1+U-U4 

i.u:)_o 03 

Q Q+0.2 

Xled (d.o.f.) 

0.999 (259) 

1.040(436) 


The parameters not indicated in this Table have been fixed to the 
best-fit values of obs.B given in the first part of Table 2 
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2.2. Phase resolved spectroscopy 

After correcting the time of arrival of the source events to the 
Solar System barycenter, we carried out a standard timing anal¬ 
ysis based on folding and phase fitting. The source period was 
clearly detected in all the observations, with the values given in 
Tabled Note that the pulse period reported in Mereghetti et al. 
(2004) for obs. B is wrong. As a consequence, also the corre¬ 
sponding pulsed fraction (computed by fitting a sinusoid to the 
background subtracted light curve) is slightly larger than that 
previously reported. The folded light curves of the three obser¬ 
vations are shown in Fig.2 (left panel). They are characterized 
by a single broad pulse of nearly sinusoidal shape and different 
levels of modulation in the three observations. We have com¬ 
puted the pulsed fraction by fitting the background subtracted 
pulse profiles with a constant plus a sinusoid. The pulsed frac¬ 
tion of each observation, defined as the ratio between the semi¬ 
amplitude of the sinusoid and the constant, is plotted as a func¬ 
tion of the source flux in the right panel of Fig. 2. 

Phase resolved spectroscopy was performed first by divid¬ 
ing the data of obs. B in 8 phase intervals as indicated in Fig. 
2. For each interval a spectrum was extracted and analyzed as 
described above for the phase averaged spectra. Again, to study 
the phase dependent variations, we adopted for simplicity the 
power law plus blackbody model. The spectra of the differ¬ 
ent intervals were first fitted keeping the parameters fixed at 
the best-fit values of the phase averaged spectrum, and allow¬ 
ing only the overall normalization factor to vary. The corre¬ 
sponding residuals and normalization factors, shown in Fig|3 
clearly indicate that the spectrum is softer at pulse minimum 
and harder at maximum. To model these spectral variations, 
as a first step we allowed the normalizations of the two spec¬ 
tral components to vary independently. Then, also the black- 
body temperature and power law photon index were left free 
to vary. The resulting best-fit values are plotted in Fig. 4. As it 
was found in the comparison of the phase averaged spectra, the 
spectral variations can be modelled equally well by a change in 
the blackbody or in the power law component. What is interest¬ 
ing is that, in any case, a change in the relative normalization 
of the two components is required. 

The same kind of analysis applied to the other two obser¬ 
vations^ led to similar results. An example of this, for the case 
of fixed kTss (fourth column of Fig. 4), is shown in Fig. 5. 

2.3. Search for spectral features 

We searched for spectral features in the EPIC spectra by first 
considering lines with a Gaussian profile and width narrower 
than the instrumental energy resolution (i.e. fixing cr=0) and 
then extending the search to broader lines (cr=100 eV). No 
spectral features were significantly detected. Since some de¬ 
viations from the best fit models, which are most likely to be 
attributed to small calibration uncertainties were present in the 
residual, we adopted a conservative confidence level of 4a for 


due to the small number of counts we used only three phase inter¬ 
vals for obs. A 


the upper limits. These are given for the equivalent widths of 
either emission or absorption lines in Table 4. 

For obs. B, thanks to the long exposure time and high 
source flux, also the RGS instrument could be used for the 
spectral analysis. The RGS covers only the 0.35-2.5 keV range 
and has a smaller effective area than EPIC, but thanks to its 
higher energy resolution (~5 eV at 1 keV compared to 80 eV 
of EPIC) it is very efficient in the search for narrow lines in 
the soft X-ray range. After a standard extraction of the RGS 
source and background spectra and the generation of the cor¬ 
responding response files, the first order spectra^ of the RGS 1 
and RGS2 instruments were rebinned in order to have at least 
50 counts per bin. The spectra were then jointly fitted using 
the blackbody plus power law best fit parameters obtained with 
EPIC, leaving only a normalization factor for each RGS unit as 
free parameter. The values of these factors are 0.99±0.04 for 
RGSl and 0.90±0.03 for RGS2, which are fully compatible 
with the known cross-calibration uncertainties. However, this 
fit is not formally acceptable (Xred =1-22/127) and some ex¬ 
cess in the residuals is visible around 1.65 keV (see Pig|6}. In 
fact, a significantly better fit (Xred =1-05/125) is obtained with 
the addition of a narrow emission line centered at 1.653)')Q gg 5 
keV and with equivalent width 21 ±7 eV (90% confidence 
level). However the reality of this line is not confirmed by the 
higher statistics EPIC data; the 4 a upper limit on the equiva¬ 
lent width of a narrow emission line at 1.653 keV in the EPIC 
spectrum of obs. B is 6.7 eV. 

Table 4. Upper limits (at 4 a) on spectral features in the EPIC 
spectra of IE 1048.1-5937 from obs. B. To derive these values 
we used narrow (cr = 0 eV) and broad (a = 100 eV) Gaussian 
lines in emission or absorption. 


Energy range 

a (eV) 

Equivalent width 

1-2 keV 

0 

< 10 eV 


100 

< 20 eV 

2-4 keV 

0 

< 20 eV 


100 

< 30 eV 

4-5 keV 

0 

< 30 eV 


100 

<55eV 

5-6 keV 

0 

< 60 eV 


100 

< 60 eV 

6-7 keV 

0 

< 90 eV 


100 

< 150 eV 

7-8 keV 

0 

< 150 eV 


100 

< 220 eV 


The search for spectral features was also extended to the 
EPIC phase resolved spectra (the smaller number of counts 
does not allow to use also the RGS for phase-resolved spec¬ 
troscopy). No narrow features were detected in the 8 phase 
resolved EPIC spectra described in Section 2.2. The 4a up¬ 
per limits for emission or absorption lines are typically a factor 
~2-5 higher than those reported in Table 4, depending on the 
pulse phase interval and energy range. 

^ the higher order spectra do not have enough counts for a mean¬ 
ingful analysis. 
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Fig. 2. Left panel: Background subtracted pulse profiles of IE 1048.1-5937 in the 0.6-10 keV energy range obtained with the 
MOSl (the only EPIC camera operated in the same mode for all observations). The phase intervals used for the pulse phase 
spectroscopy of obs. B are indicated. Right panel: Pulsed fraction as function of the average count rate. 


3. Discussion 

The target of opportunity observation of July 2004 (obs. C) was 
motivated by the detection of one burst from the direction of 
IE 1048.1-5937 with the RXTE satellite (Kaspi et al. 2004) 
and carried out nine days after this event. No bursts were seen 
during the XMM-Newton observation, nor any other sign in¬ 
dicating a particularly increased level of activity: the flux of 
IE 1048.1-5937, 6.2x10-12 erg cm-^ s"!, was still ~25% 
higher than the typical values observed before the 2001 out¬ 
burst (see Eig.2 in Mereghetti et al. 2004 and Eig.l in Gavriil & 
Kaspi 2004), but not as high as that measured in obs. B. We find 
statistically significant differences in the spectra of the three 
XMM-Newton observations, however such differences are not 
related in a monotonic way with the source luminosity: when 
the flux was at the highest level (obs. B) the spectral hardness 
was intermediate (see Eig.l). 

A coherent pattern is instead present in the relation between 
flux and pulsed fraction: the latter decreases when the source 
brightens, as was first reported by Mereghetti et al. (2004) 
based only on two observations. This is confirmed by obs. C 
as shown in the right panel of Eig. 2, where, to avoid any cross 
calibration uncertainty, we directly compare the flux and pulsed 
fraction (Pp) measured with the same detector in the three ob¬ 
servations. The existence of an empirical anti-correlation be¬ 
tween luminosity and pulsed fraction, independent on its theo¬ 
retical explanation, must be taken into account in the interpre¬ 
tation of the RXTE results for this source. The luminosity and 
fluence obtained with this satellite, which can only measure the 
pulsed component of the flux, are clearly underestimated if the 


decrease in pulsed fraction during the outbursts is not taken 
into account. We estimate that total energy release of the flares 
peaking in November 2001 and June 2002 is at least the double 
of the value (2 and 20 x 10"^° ergs respectively) derived assum¬ 
ing a constant Pf=. 94 by Gavriil & Kaspi (2004). 

Our results for the phase averaged spectroscopy are sub¬ 
stantially in agreement with previous findings for this and other 
AXPs, except that the high statistics of obs. B allowed us to re¬ 
ject the fit with two blackbody components. We also explored 
spectral fits using a simple Comptonization model. The fit with 
the Comptonization model alone is not satisfactory. The addi¬ 
tion, however, of a second, softer blackbody component, re¬ 
sulted in an acceptable fit, both for the phase-averaged, and for 
the phase-resolved spectra. Quite interestingly the estimate of 
the radiating area for the cooler component is in this case com¬ 
patible with the entire star surface when the distance uncer¬ 
tainty is taken into account. Spectra at different phases can be 
fitted keeping this emitting area fixed at the value derived from 
the phase-averaged analysis and this results only in a modest 
variation of the cool blackbody temperature, « 10%, and a 
modulation of the hot blackbody emitting area. The smaller 
area associated with the hotter blackbody emission may be sug¬ 
gestive of a scenario in which a magnetically active region pro¬ 
duces a hotter region contributing a substantial part of the lu¬ 
minosity. The accelerated high-energy particles, besides heat¬ 
ing this region may upscatter the soft photons giving rise to the 
Comptonized spectrum. 

The implementation of a detailed Comptonization model 
to fit the X-ray spectrum is beyond the scope of our anal¬ 
ysis. Therefore, we have used the approximated model de- 
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Fig. 3. Phase resolved spectroscopy of obs.B. The figures show the ratio between the phase resolved spectra and the best-fit 
power-law plus blackbody model rescaled by a factor F (indicated in each panel). The data (black PN, red MOSl, green MOS2) 
have been graphically rebinned to emphasize the spectral trend. 


scribed in Section 2.1, which is based on a constant scatter¬ 
ing cross-section. The scattering cross-section in a strongly 
magnetized medium is anisotropic and quite different for or¬ 
dinary (O) and extraordinary (X) photons (e.g. Ventura 1979; 
Gonthier et al. 2000). In our case the typical incident pho¬ 
ton energy is « 0.5 keV, well below the resonance at the 
electron cyclotron energy Eb ~ 11.6(B/10^^ G) keV. In the 
low-frequency limit (E <IC Eb), the cross section for ordinary 
photons is ao ctt sin^ 6, where ut is the Thomson cross- 
section. Due to relativistic beaming, scattering occurs at an in¬ 
cident angle 9 ^ 1 /( 27 ) in particle frame and therefore 
ao ~ crT/(47^). For relativistic electrons of a given energy, 
the cross section for ordinary photons is therefore constant. 
The cross section for extraordinary photons, instead, peaks at 
9 ^ (E/Eb)^^^, where ctx ~ crT(E/EB). It means that, tak¬ 
ing E/Eb « 10“^, it dominates over ctq for 7 ^ 15. In this 
case the use of a constant cross-section is possible because we 
have checked numerically that its dependence on the photon 
spectrum has a negligible effect on the Comptonized spectrum. 

The present model does not allow to univocally determine 
the properties of the Comptonizing medium, since the fitted pa¬ 
rameter a depends on both the scattering depth and the particle 
energy. However, with the presently derived values of the pho¬ 
ton index, a ~ 5, once the Comptonizing particles are assumed 
to be relativistic electrons, the medium is optically thin and the 
optical depth is strongly anticorrelated to the particle energy 
(tS ~ 7~®)- For example, if 7 ^ 5, the magnetic depth is 

~ 10“®, to which corresponds a non-magnetic Thomson 


depth ^ 10“^, while if 7 ~ 50, drops to 10“corre¬ 
sponding to a non-magnetic Thomson depth of ~ 10“^^. 

Thanks to the high statistics of our observations we could 
set strong upper limits on the presence of spectral lines (in view 
of the EPIC results we do not consider the possible feature at 
1.65 keV seen in the RGS as really present in the source). If in¬ 
terpreted in terms of magneto-dipolar braking, the spin-down 
rate measured with RXTE at the time of obs. B (Gavriil & 
Kaspi 2004) yields a dipolar field Ri 8 x 10^^ G and at these 
field strenghts cyclotron, free-bound and free-free features are 
expected to lie in the observed energy band. On the other hand, 
present estimates of their detectability in the presence of su¬ 
perstrong magnetic fields still suffer from a number of approx¬ 
imations (Zane et al. 2002, Ho et al. 2003). Several effects, 
e.g. fast rotation and magnetic smearing, have been suggested 
to explain the absence of spectral lines in other sources (e.g. 
Braje & Romani 2002). Although the former does not apply to 
the present source, both cyclotron and atomic lines are sensitive 
to the local field strength and the contributions of surface ele¬ 
ments with different values of B may under certain conditions 
wash out spectral features. Ho & Lai (2004) also suggested 
that at super-strong field strengths {B ^ 10^^ G) vacuum po¬ 
larization may efficiently suppress spectral lines, but current 
treatment of this QED effect is too crude to make a definite 
statement. The absence of lines may also suggest either that the 
magnetic field at the neutron star surface is larger than the large 
scale dipolar component (i.e. that i? > 15. x 10^"^ G) or, alter¬ 
natively, that the value inferred by the spin-down rate is over- 
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Fig. 4. Best-fit parameters for the phase resolved spectra of obs. B. From top to bottom; the observed 2-10 keV flux (in units of 
10“^^ erg cm“^ s“^), the ratio of the 0.6-10 keV (absorbed) flux in the blackbody and power law components, the power law 
photon index, the blackbody temperature and the reduced The diamonds indicate that the parameter is fixed to the best-fit 
value of the phase averaged spectrum. All the error bars represent Icr errors. 


estimated (i.e. that B < 10^"^ G). As discussed by Thompson 
et al. (2002), this value must be regarded as an upper limit: a 
twisted magnetosphere may lead to a reduction up to an order 
of magnitude in the inferred polar value of the magnetic field, 
while the estimate of the spin down age is unaffected. 

Phase resolved spectroscopy has now been carried out for 
most AXPs. The earlier observations obtained with the ASCA 
and BeppoSAX satellites gave some indications for phase- 
dependent spectral variations in 4U 0142-1-614 (Israel et al. 
1999), IE 22594-586 (Corbet et al. 1995, Parmar et al. 1998), 
IE 1048.1-5937 (Oosterbroek et al. 1998), and a more sig¬ 
nificant evidence for IRXS J1708-4009 (Israel et al. 2001). 
These indications were later confirmed by observations with 
higher statistics carried out with XMM-Newton (Tiengo et al. 
2002, Woods et al. 2004, Gohler et al. 2004) and Chandra 
(Morii et al. 2003). In most of these studies the blackbody plus 
power-law was adopted and the spectral variations described 
as changes in one or more of the model parameters. The re¬ 
sults reported here give strong evidence that the spectrum of 
IE 1048.1-5937 is harder at pulse maximum. However the in¬ 
terpretation of this result is not straightforward since the val¬ 
ues of the spectral parameters, e.g. blackbody temperature or 
photon index, are linked in a complex way to the overall spec¬ 
tral shape/hardness. As we have shown in the previous section, 
these variations can be modelled in different ways. Eor exam¬ 
ple (see EigH if the photon index and blackbody temperature 
are fixed, the ratio Fbb/Fpl is proportional to the total inten¬ 
sity. In turn, this means that the blackbody contribution at the 
pulse maximum is more important than that at the pulse mini¬ 


mum, i.e. that the blackbody component has a pulsed fraction 
larger than the power law. This is what is expected if the ther¬ 
mal component originates from the star surface and the non- 
thermal one from a more extended magnetosphere, and may 
lead to suggest that the observed spectral changes are simply 
due to the partial occultation of an the emitting regions as the 
star rotates, with the effect being more prominent for the sur¬ 
face emission than for the magnetosphere. On the other hand, 
better fits are obtained if also the photon index or the tempera¬ 
ture are free to vary (third and fourth columns of Eig. 4). In this 
respect, probably the most interesting result is that we found 
impossible to reproduce the phase-dependent variations keep¬ 
ing one of the two components completely hxed (i.e. shape and 
normalization). This fact indicates a coupling between the two 
components, or alternatively, it strengthens a scenario in which 
the emission below ~10 keV is dominated by a single physical 
mechanism with a spectrum more complex than a blackbody or 
a power-law. 

4. Summary 

The main conclusions derived from the analysis and compar¬ 
ison of the three XMM-Newton observations of IE 1048.1- 
5937 reported here are the following: 

- long term flux variations of a factor ~ 2 are accompanied 
by only minor, though statistically significant, changes in 
the source spectral properties 

- the pulsed fraction is anti-correlated with the source lumi¬ 
nosity, decreasing from ~90% to ~55% in correspondence 
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Fig. 5. Best-fit parameters for the phase resolved spectra. The horizontal lines indicate the values for the phase averaged spectra. 
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of a twofold increase in flux. This should be properly taken 
into account when the source energetics are inferred by 
measurements of the pulsed flux 

- spectral variations are clearly present as a function of the 
pulse phase, with the hardest spectrum at pulse maximum; 
their interpretation in terms of simple changes in two com¬ 
ponent models is not unique 

- there is no evidence for absorption and/or emission lines in 
the spectrum with equivalent width larger than 10-220 eV, 
depending on energy and width 

- the spectrum is consistent with a two component model 
consisting of a blackbody emission from an emitting area 
compatible with the whole neutron star surface and a hotter 
Comptonized blackbody, possibly coming from a relatively 
large region of the star surface heated by non-thermal mag- 
netospheric particles also responsible for the scattering. 
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Fig. 6. First order RGSl (black) and RGS2 (red) spectra of IE 1048.1-5937 during obs. B. The sharp drops in the data and model 
are due to the gaps between chips. No RGSl data are present in the ~0.9-1.2keV range because of the failure of one CCD. The 
lower panel shows the residuals (in standard deviation units) with respect to the EPIC best-fit power law plus blackbody model. 








































































































































































































